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The excitonic absorption spectrum for a disordered quantum well in presence of exciton-acoustic 
phonon interaction is treated beyond the Markov approximation. Realistic disorder exciton states 
are taken from a microscopic simulation, and the deformation potential interaction is implemented. 
The exciton Green's function is solved with a self energy in second order Born approximation. 
The calculated spectra differ from a superposition of Lorentzian lineshapes by enhanced inter-peak 
absorption. This is a manifestation of pure dephasing which should be possible to measure in 
near-field experiments. 

PACS numbers: 78.20.-e, 63.20.Ls, 73.43.Cd 
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INTRODUCTION 



Unique information on semiconductor nanostructures 
even beyond the diffraction limit can be obtained using 
near-field spectroscopy. This opens the way to detect 
fine spectroscopic and structural details 1 . Usual selec- 
tion rules for optical transitions are broken, which al- 
lows one to observe also states which are dark in the 
far-field 2 . The present record experiment is a spatial res- 
olution of 30 nm which was achieved using a tapered 
optical fiber and an InAs self-assembled quantum dot 
(QD) sample 3 . Even with a less demanding setup us- 
ing a resolution comparable to the wavelength of light, 
individual optical transitions can be seen due to their 
spatial and spectral position. Absorption spectra from 
quantum wells (QW) with interface fluctuations could 
be measured in transmission geometry through a metal 
aperture of 400 nm on a GaAs/AlGaAs sample whose 
substrate was removed 4 . This gives access to single exci- 
ton state absorption lineshapes, with the chance to detect 
deviations from a Lorentzian profile. Similar lineshapes 
are seen in photoluminescence (PL). Among others, Be- 
sombes et al. 5 found that the PL from a single CdTe QD 
exhibits broad bands on the sides of the Lorentzian peak. 
This has been also seen in Fourier-transformed four wave 
mixing experiments by Borri et al. 6,7 , who systematically 
studied the importance of the broad band with respect 
to the (Lorentzian broadened) zero phonon line (ZPL) as 
a function of confinement energy in InAs/InGaAs QDs. 

However, using the standard Fermi's golden rule for 
the phonon scattering, only a Lorentzian broadened ZPL 
is expected for the absorption spectrum of a single exci- 
ton state. Phenomena which go beyond the framework 
of Fermi's golden rule (such as the phonon broad band) 
are often called quantum kinetic or non-Markovian (i.e. 
memory) effects. The broad band can be understood in 
terms of a non-perturbative coupling between electronic 
and phonon degrees of freedom. This goes back to the 



Huang-Rhys theory of -F-centers 8 , which considers intra- 
state (diagonal) coupling to phonons only. This problem 
can be solved exactly with the cumulant expansion (also 
known as Independent Boson Model 9 ). Since off-diagonal 
level coupling is absent, the ZPL remains unbroadened. 
Within this model, different exciton-phonon couplings 
such as optical polar, acoustic piezoelectric, and defor- 
mation potential have been studied 10 . This model was 
recently extended to include off-diagonal coupling, suc- 
cessfully explaining why the ZPL get broadened even if 
real phonon-assisted transitions to other levels are not 
possible (pure dephasing) 11,12 . Apart from the cumulant 
expansion, approaches based on various dynamical ap- 
proximations have been used, too. For instance, a 4th 
order selfconsistent summation of diagrams allowed one 
to compare QW linewidths extracted from photon echo 
measurements 13,14 : It was found that at elevated temper- 
atures the linewidth broadening is mainly due to pure 
dephasing. However, a simple quantum disk model for 
computing the exciton wave functions was used in this 
case. The correlation expansion for the carrier-phonon 
interaction which was used in Ref.15 to model the cou- 
pling to both optical and acoustical phonons, is exact up 
to the 2nd order. However, neither Coulomb or disorder 
effects were considered there. 

A consistent theoretical treatment which describes 
non-Markovian effects for realistic, disordered excitonic 
resonances from semiconductor QWs is, to the best of 
our knowledge, still missing. In this paper we face this 
problem accounting for both diagonal and off-diagonal 
coupling to acoustic phonons via deformation potential 
interaction. The excitonic self energy is taken in sec- 
ond Born (2B) approximation and integrated numeri- 
cally. This approach reproduces both the real scattering 
processes leading to the major ZPL broadening and the 
level-diagonal pure dephasing responsible for the broad 
band. The resulting (near-field) absorption spectrum ex- 
hibits apart from the Lorentzian-broadened individual 
lines a superposition of broad bands which results in an 
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enhanced inter-peak absorption. We suggest that this ef- 
fect of non-Markovian dynamics can be demonstrated in 
near-field measurements. 

The paper is organized as follows: The theory with 
a presentation of different approximations for the self 
energy in Sect. (II) is followed by numerical results in 
Sect. (Ill), and the conclusions are given in Sect. (IV). 
Finally, two appendixes contain more technical details. 



II. 



THEORY 



To model a realistic QW with disorder it is convenient 
to start with disorder eigenstates of the exciton. If the 
disorder strength is much smaller than both the level 
distance of electron and hole confinement states and the 
exciton binding energy, the total exciton wave function 
(in envelope function and effective mass approximation) 
can be factorized as 16 

* a (r e ,r h ) = u e (z e ) U h (Zh) <t>u{p e - Ph)^a(R) , (1) 

where confinement u a (z a ) in growth direction z, in-plane 
exciton relative motion <pi s (p), and in-plane center-of- 
mass (COM) wavefunction ip a (H) appear. Since the first 
two functions are independent of disorder, they have to 
be computed just once for the disorder-free QW. Then, 
it suffices to solve a single particle equation for the COM 
exciton states in two dimensions, 



V>a(R) = € a i> a {K). (2) 



The potential V(R) accounts for interface disorder and 
is spatially correlated on the scale of the exciton relative 
motion, 



V(R) 



a—e,h 



(3) 

with mass ratios r\ e = M/rrih and rjh — M/m e . The well 
width fluctuations AL Z (R) transfer into potential fluc- 
tuations via the derivative of the ideal electron and hole 
confinement energy dE a /dL z 17 . As zero of energy, we 
take the ideal Is exciton transition energy of the lowest 
sublevel transition hhl-el in a GaAs/Al .3Gao.7As QW 
of given width L z . 

Recently, Grochol et al. 18 have critically questioned 
the factorization Ansatz Eq.(l). Calculations with the 
full in-plane electron-hole wave function gave systemati- 
cally lower eigenergies, while the wave functions showed 
only minor modifications compared to the factorization, 
at least in the low-energy part of the spectrum. Since 
in the present paper, we arc interested in the low-energy 
tail of the absorption spectrum with well localized states, 
the factorization can be safely used. 

The interaction with acoustic phonons is treated here 
as the main source of scattering and dephasing. Optical 



phonon scattering is not important here since, at temper- 
atures below room temperature, optical phonon absorp- 
tion by the exciton states is unlikely. Further, piezoelec- 
tric scattering is neglected since the exciton wave func- 
tions arc typically large in momentum space, leading 
to small matrix elements for this kind of coupling 19,20 . 
Thus, only deformation potential interaction with longi- 
tudinal acoustic phonons is relevant. 

This model leads to the following Hamilton opera- 
tor for excitons (creation operators and acoustic 
phonons (creation operators a q with three-dimensional 
momentum q) with (bulk like) dispersion u> q = us |q|) : 



U = ^2e a BlB a + ^ftw q a q a q 

a q 

a/3q 



(4) 



The deformation potential matrix elements t^a are dis- 
cussed in App. (A). We have neglected the spin degree of 
freedom which splits each localized exciton state into a 
bright doublet via the exchange interaction. This split- 
ting is driven by COM anisotropy and has typical values 
of a few tens of /zeV 17 . It can be disregarded in calculat- 
ing absorption spectra with not too high resolution. 

Within linear response theory, the optical properties of 
the system are given by the one-exciton Green's function 
(exciton propagator) 



G af3 (t-t') = -i(TB a (t)BUt>)) 



(5) 



with time ordered finite-temperature expectation value. 
In frequency space, the Green's function obeys the Dyson 
equation (in what follows h = 1 is used, and the frequency 
u> is understood to carry a small negative imaginary part, 
w-iO) 



[(W - £ a )5 a ri - EchjM] G n p(uj) = S a p . 



(6) 



introducing the exciton self energy matrix H arj {u)). The 
absorption spectrum is proportional to 



a(u) = Im^m* G a p(w) mp . 

a/3 



(7) 



The optical interband matrix elements m a arc defined in 
App. (A). 

For the self energy, we take the lowest order diagram 
for the exciton-phonon interaction shown in Fig. 1. It 
is called 2nd order Born (2B) approximation since the 
interaction appears just twice: 



q 



1 + n(ui q ) 



n(Wq) 



LU - 6p - Wq U) - Ep + Wq 



(8) 

The phonons arc assumed to be in equilibrium at the lat- 
tice temperature T, n(w q ) = l/(exp(u> q //cBT) — 1). The 



two terms in Eq. (8) are related to phonon emission and 
absorption. Note that the internal line in Fig. 1 is under- 
stood to represent the bare exciton propagator. Taking 
instead a dressed propagator with shifted quasiparticle 
energy would be an improvement. Our calculations have 
shown that this (acoustic) polaron shift for the localized 
exciton states is of the order of a few /ieV only, and we 
have chosen to neglect it everywhere. 

An alternative derivation would start with the equa- 
tions of motion for the exciton propagator which cou- 
ples to mixed exciton-phonon operator expectation val- 
ues (phonon-assisted density matrix). The next equa- 
tion in the hierarchy is then decoupled in a way that 
two phonon operators can be combined into the phonon 
occupation function. Eliminating the phonon assisted 
density matrix from the equations leads to a differen- 
tial equation of the exciton propagator with a scattering 
being essentially a time integration over the past. This 
kind of memory in the scattering is fully equivalent to 
the frequency dependence of the self energy in the di- 
agrammatic approach. Further, the non-diagonality of 
the self energy S a ^(uj) leads to cross correlations among 
individual excitons. 

In the Markov approximation these correlations are 
neglected, but more importantly the frequency depen- 
dence of the self energy is dropped completely. In the 
time frame, this means to extract the exciton propaga- 
tor from the scattering integral at the latest time, thus 
skipping any memory effect. Consequently, the tempo- 
ral dynamics of the propagator reduces to an exponen- 
tial decay with a constant damping j^f . In this paper 
we go beyond this level of approximation, aiming at a 
full non-Markovian description which includes the non- 
diagonality and the memory (frequency dependence) in 
the self energy. For the sake of comparison, we discuss in 
the following subsections the solution of Eq. (6) at various 
levels of sophistication. 

The radiative contribution to the self energy (radiative 
damping) will be neglected since it is tiny in thin QWs 
which have strongly localized excitons 17 . 
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FIG. 1: Second order self energy diagram. The arrow denotes 
the exciton propagator, while the wiggly line stands for the 
phonon propagator. Each dot carries a coupling element t. 



decides which exciton states are coupled to the phonon 
bath. It starts at least with a power of E 3 at small 
energies, as seen in App. (A). Since the exciton wave 
functions are real valued, the following symmetries hold: 
J^{E) = J^(E) and J&,(25) = Jfo{E). Thus, if we 
consider a system of A exciton states, (A 3 + N)/2 com- 
ponents of the coupling function have to be computed. 



B. Markov approximation 

The Markov approximation for the exciton-phonon dy- 
namics is defined within our formalism by keeping only 
diagonal terms of the self energy and evaluating them at 
the exciton resonance (on-shell), 

£ a „(w) - d aT ,X™ , = H aa {w = e a -iQ). (11) 

Apart from the (small) polaron shift originating from the 
real part, this level of approximation implies a constant 
damping 

7^/2 = ImE^ =nJ2 n(ep - e a ) jfje, - e a ) . (12) 





A. Full non-Markovian level 

The self energy Eq. (8) can be written in a more com- 
pact form by extending the energy integration to positive 
and negative values, 

r J p (E) 

Sa„M=£ dEn(E) a f> (9) 

J U! - 6/3 + E 

which accounts for both phonon absorption and emission. 
The coupling function 

J^(E) ^sgn(E)J2t^5(\E\ - Wq ). (10) 
q 

can be viewed upon as a phonon density of states 
weighted by coupling matrix elements. Its spectral range 



(The factor 2 appears because 7^ is the decay rate of the 
intensity, which is proportional to the squared exciton 
propagator.) Thus, Eq. (6) reduces to a simple algebraic 
equation, and the absorption comes out as 

which is a Lorentzian lincshape for each line. No broad 
bands are obtained within the Markov approximation. 
The rate 7^ is identical to the result of Fermi's golden 
rule. We notice that within the Markov approximation, 
only (A 2 + A) /2 components of the coupling function 
have to be computed. 
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C. Single state limit 

As already mentioned, it is not possible to account for 
broad bands within the Markov approximation. In the 
single-state limit (uncoupled exciton states), the Markov 
rates Eq. (12) vanish since the coupling function Jg a (E) 
tends strongly to zero at zero argument E = e a — e a . 
Thus, in the absence of radiative damping, the absorption 
consists of delta lines. However, if we do not perform the 
Markov approximation, a composite absorption lineshape 
is obtained due to the frequency dependence of the self 
energy T, aa (u). A broad band adds to the unbroadened 
zero phonon line. Its weight with respect to the total 
absorption area is Z a = 1/(1 + S a ), with the Huang- 
Rhys factor S a stemming from the self energy derivative 

A« = j = / dE n(E) — — — > . 

(14) 

This result for Z a agrees up to first order in S a with the 
exact result from the independent Boson model, Z a = 
exp(— S a ). However, the single-state solution neglects 
completely the inter-state scattering which leads to the 
Lorentzian broadening of the lines. 

D. Semi-Markov approximation 



100 




x [nm] 

FIG. 2: (Color online) Probability amplitude of six exciton 
COM wave functions in a 5 nm wide GaAs/AlGaAs QW. The 
lines are drawn at 50, 30, and 10% of the peak probability 
of each wave function. Notice that state 1 and state 6 do 
overlap. 



In the full non-Markovian solution, both the 
Lorentzian broadening of the ZPL and the broad bands 
are obtained using the full self energy Eq. (9). However, 
the numerical cost of 0(N 3 ) coupling functions makes 
this treatment unfeasible for the large number N of ex- 
citons states covered by the focus of a typical near-field 
experiment. Thus, either one has to select fewer states, 
or one has to resort to another approximation which still 
preserves the main features of the absorption spectrum. 

Aiming at this goal, we introduce the "Semi-Markov 
approximation" which neglects the off diagonal terms in 
the self energy and keeps the frequency dependence only 
in the fully state diagonal terms: 



n(E)JS a (E) 



dE- 



E 



uj — e a + E 

n{E)jP a {E) 



dE- 



iO 



(15) 



A linear expansion of Sf M (w) around the exciton pole 
e Q shows (cf. App. B) that the ZPL has a reduced width 



7f M = 7*7(1 + S a ) 



(16) 



The semi-Markov approximation is acceptable as long as 
the off diagonal elements of the self energy are small with 
respect to the differences between diagonal elements of 



the matrix problem Eq. (6). Since these are dominated 
by the exciton eigenergies, the semi-Markov approxima- 
tion is expected to deviate markedly from the full non- 
Markovian treatment only in spectral regions with a high 
density of overlapping states. Although much less de- 
manding in computer power (only (TV 2 + N) /2 coupling 
functions needed - as in the Markov approximation) , both 
broad bands and broadening of the ZPLs are obtained in 
the semi-Markov approximation. 



TABLE I: Material constants for bulk GaAs and other pa- 
rameters used in the numerical simulations. 



longitudinal sound velocity 


us 


5.33 


nm/ps ° 


bulk mass density 


PM 


5.37 


g/cm 3 " 


def. pot. electrons 


A 


-7000 


meV a 


def. pot. holes 


D v 


+3500 


meV a 


exciton mass 


M 


0.300 


m » 


in-plane electron mass 


m e 


0.067 


m 


in-plane hole mass 


rn h 


0.233 


m " 


QW width 


L z 


5 


nm 


disorder strength 


a 


4.75 


meV 


correlation length 


(IB 


9.9 


nm 


grid size 




1.65 


nm 



"After Ref. 21 
6 After Ref. 22 
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FIG. 3: (Color online) Selected components of the coupling 
function C% V (E) = Ji v (E)/E 3 for some of the states of Fie. 2. 
The off-diagonal components have been magnified by the fac- 
tors given in brackets. 



III. RESULTS 

We present now numerical solutions of the the- 
ory developed in Sect. (II). The Schrddinger equation 
Eq. (2) has been discretized on a square grid with step 
size A x . An uncorrelated potential corresponding to 
(dE a /dL z )AL z (R,) in Eq. (3) has been generated and 
convoluted with the relative exciton wave function, us- 
ing parameters of a 5 nm wide GaAs/Alo.3Gao.7As QW. 
The resulting disorder potential V(R) is computed on 
the grid and has a correlation length close to the exciton 
Bohr radius as of this specific QW. The variance a of this 
correlated potential has been adjusted to fit experimental 
results from the speckle analysis 23 . The diagonalization 
of Eq. (2) on a simulation mesh of 128 x 128 grid points 
with periodic boundary conditions was achieved via an 
ARPACK-based package. All used material constants 
and simulation parameters are listed in Tab. I. In partic- 
ular the in-plane hole mass is taken from an excitonic k -p 
calculation by Siarkos et al. 22 which takes into account 
the heavy-light hole mixing. An effective hole mass of 
rrih = 0.233 mo could be extracted which describes sat- 
isfactorily the exciton COM motion for a wide range of 
quantum well thicknesses. We have used this hole mass 
for the effective-mass calculation in the present paper. 

In Fig. 2 we display the amplitudes |^; Q (R)| 2 of 
six computed states, whose eigenergies are contiguous. 
These states were selected by purpose of exhibiting vari- 
ous combinations of oscillator strength and spatial over- 
lap to other states. In particular both quantities are large 
for ipi(R) and both are small for ^(R). State ^s(R) 
has large oscillator strength and small overlap, while for 
-06 (R) the contrary holds. Here, we are restricted to six 
states only since we are going to compare the computa- 
tionally most demanding spectra using the full self energy 
matrix Eq. (9) with other approximations. 



FIG. 4: (Color online) Absorption lineshape for a single state 
at T = 40 K (upper curves, red) and T = 10 K (lower curves, 
blue). Results from the 2nd Born approximation (solid lines) 
are compared with the exact solution (dashed) for state a = 1 
(left panel) and a = 3 (right panel) of Fig. 2. The localization 
lengths are Ai = 16.4 nm and A3 = 23.3 nm. The arrows 
indicate the ZPLs, which correspond also to the Markov limit 
of the absorption spectra. 

The coupling function J^ n (E) is formed using the 
eigenstates ^> a (R) and represents the most time- 
consuming part of the numerical task. The fully diagonal 
components J^ a (E) are responsible for pure dephasing in 
the single-state limit. These components, divided by E 3 , 
are approximately Gaussians with FWHM inversely pro- 
portional to the localization length A a of state a (full 
curves in Fig. 3 for a = 1 and a = 5). Their local- 
ization lengths (calculated from the inverse participation 
ratio) are 16.4 and 36.9 nm, respectively. Partially off- 
diagonal terms J^ a (E) with (3 ^ a appear in the Marko- 
vian scattering rates 7^ (Eq. (12)) and are sensitive to 
the spatial overlap between ip a (R) and ^g(R). In this 
sense, V'i(R) overlaps with ipQ (R) about 100 times more 
effectively than with i/> 5 (R), according to Fig. 3. Finally, 
fully off-diagonal terms J^ n (E) with a, f3, rj all different 
account for deviations between the full solution with self 
energy Eq. (9) and the semi-Markov level with self energy 
Eq. (15). All nondiagonal terms of the coupling function 
vanish at E — > with power E 4 due to the orthogonality 
of the exciton wave functions. 

Having computed the coupling functions, the self en- 
ergy at the desired level of approximation (single state, 
Markov, semi-Markov, or full non-Markov ) can be easily 
generated. Using this input, the Green's function and the 
absorption spectrum are obtained. First we have checked 
the quality of our 2B treatment against the exact solu- 
tion (Independent Boson Model) which is available in 
the single-state case only. In Fig. 4 we show the absorp- 
tion spectra for state a = 1 and a = 3 of Fig. 2. The 
broad band stems from phonon absorption (left side of 
the ZPL) and phonon emission (right side). Therefore, 
the asymmetry evolves into a more symmetric shape as 
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FIG. 5: (Color online) Absorption spectrum calculated in 
the Markov approximation (black dashes), with semi-Markov 
quality (blue dots), and with full non-Markovian quality (red 
full line) at lattice temperatures of T = 10 K (top panel) and 
T — 40 K (bottom panel). Notice the good agreement of 
semi-Markov and full non-Markovian results. No focus func- 
tion is included. The triangles indicate the positions of the 
bare exciton eigenergies. 

temperature is raised. There is very good agreement be- 
tween 2B and exact results at T = 10 K. However, at 
higher temperatures the broad band is underestimated 
in the self energy approach. 

In Fig. 5 we display the absorption spectrum for the 
system of six exciton sates. The Markovian spectrum 
(dashed) is simply a superposition of Lorentzian lines, 
whose widths are set by Fermi's golden rule, Eq. (12). 
Since the deformation potential matrix elements t^p are 
sensitive to the spatial overlap of the COM wave func- 
tions (see App. (A)), states 2,3,5 exhibit narrower lines 
compared to states 1,4,6 which have a stronger overlap. 
The full non-Markovian spectrum is characterized by en- 
hanced absorption between the single lines, as stressed 
by the logarithmic representation of Fig. 5. Furthermore 
a shift of the peak energies is observed (polaron shift). In 
the low energy part of the spectrum (E < —3 meV) the 
energetic spacing among the states is comparable to or 
larger than the width of diagonal elements of the coupling 
function (cf. Fig. 3). Thus, sidebands are clearly noticed 
although they overlap. They are about one order of mag- 
nitude larger than the superposition of the Lorentzian 
tails in the Markovian spectrum. Conversely, the ener- 
getic distance between states 5 and 6 is much smaller 
than the width of the corresponding coupling function 
J° n (Ji). In this region the inter-peak absorption is dom- 
inated by the Lorentzian tails of the ZPLs, and the in- 
dividual broad bands cannot be easily resolved. A care- 
ful inspection reveals that the Lorentzian peak widths 
are smaller than in the Markovian case. According to 
Eq. (16) and Eq. (14), the width reduction scales with the 
Huang-Rhys factor and is more important for stronger lo- 
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FIG. 6: (Color online) Temperature dependent ZPL widths 
7« for lines 1 and 3 of the full non-Markovian (solid red lines) 
and the semi-Markov spectrum (dotted blue lines) of Fig. 4. 
The simplified expression Eq. (16) (black solid lines) and the 
Markov rates 7^ (black dashed lines) are also displayed for 
comparison. 



calized states and higher temperatures, as observed (cf. 
e.g. the peak of state 1). In Fig. 5 we also check the qual- 
ity of the semi-Markov approximation against the full 
solution. The agreement is excellent at both T = 10 K 
and T = 40 K. Thus, pure dephasing is dominated by 
state-diagonal processes. A similar conclusion has been 
obtained in Ref. 14. 

In Fig. 6 the temperature dependence of the linewidths 
of two selected states of Fig. 5 is studied. The reduced 
Markov width Eq. (16) fairly describes both full solution 
and semi-Markovian widths, as extracted from Lorentz 
fits in a narrow spectral window around the ZPLs. This 
linewidth reduction from Markov to non-Markov is a gen- 
eral feature and appears as soon as the spectral weight 
of the quasi particle is reduced. However, being a higher 
order effect in the coupling strength, one may ask how an 
inclusion of higher order terms would alter this finding. 
A consequent treatment of all diagrams up to 4th order in 
the cumulant expansion by Muljarov et al. 12 has shown 
that the reduction is indeed a pertinent feature. Its nu- 
merical value, however, might be quantitatively not cor- 
rect when restricting to the level of 2nd Born approxima- 
tion. In Fig. 6, a sublinear dependence of the linewidth 
on temperature is found for both states. This can be 
easily explained by noting that the Huang-Rhys factor 
increases nearly linear with T (cf. App. (B)). 

Having tested the features of the semi-Markov approx- 
imation in Fig. 5 and Fig. 6, we apply it to a system of 
many exciton states. To come closer to the experimental 
situation of a near-field measurement, we have included 
a Gaussian focus of 50 nm width (FWHM of field inten- 
sity) . Using the 50 lowest energy eigenstates of a disorder 
realization for a 5 nm wide QW we obtain the spectrum 
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FIG. 7: (Color online) Top panel: ZPL weights Z a = 1/(1 + 
S a ) (blue squares) and A a /(A ol +bT) ratio (black circles) for 
a set of 50 exciton states at T = 40 K. A good agreement is 
predicted from Eq. (B5) and it is found using b = 0.12 nm/K. 
The dashed line is the constant ZPL weight within the Markov 
approximation. Bottom Panel: Markovian (shadowed area) 
and semi-Markovian (blue solid line) absorption spectrum for 
the same set of states at T — 40 K using a focus of 50 nm. 



in Fig. 7 up to E = meV which is the exciton energy in 
an ideal QW of the same width. In a near field geometry, 
the visibility of a state depends not only on its oscillator 
strength, but on the overlap with the illumination spot 
as well, as seen from Eq. (A2). Thus, several states in the 
range —3 meV< E < —2 meV appear especially strong in 
the present simulation. The absorption background be- 
tween the peaks can be clearly seen on a linear scale in 
this spectral region. In the high energy tail (E « meV) 
there is practically no difference between Markov and 
semi-Markov spectra. Even on a logarithmic scale both 
spectra are nearly identical (not shown) . These observa- 
tions can be explained considering the localization length 
A a and the spectral separation of the exciton states. A a 
is increasing across the spectrum from A2 = 17.2 nm (the 
ground state is dark in the actual illumination geome- 
try) up to A50 = 54.1 nm. Via the Huang-Rhys factor, 
this leads to an increasing ZPL weight (cf. App. (B)) 
which goes from Z2 — 0.77 up to Z50 = 0.92, as seen 
in the top panel of Fig. 7. Thus, the importance of the 
non-Mar kovian broad bands (proportional to 1 — Z a ) de- 
creases with photon energy. However, the observability 
of this effect strongly depends on the spectral separation 
of the resonances, as discussed for Fig. 5. The spectrally 
close doublet at about —4.2 meV e.g. does not show a dis- 
tinct broad band, in spite of having a short localization 
length (large S a ). Thus, for observing non-Mar kovian 
effects in the absorption spectrum from QWs, it is nec- 
essary to position the near-field focus in such a way that 
states with small localization length and large spectral 
separation are excited. 



A non-Markovian computation of the absorption which 
includes all optically active exciton states is at the mo- 
ment beyond our computational possibilities. We expect 
somewhat broader ZPL peaks since more states would be 
available for phonon-mediated scattering. However, the 
broad band absorption between peaks is not expected to 
be modified as it is due to intra-state relaxation (first 
term on the rhs of Eq. (15)). Thus, a calculation of the 
absorption spectrum using all states will not reveal new 
features as compared to what we display in Fig. 7, at least 
within the semi-Markov approximation. 



IV. 



CONCLUSIONS 



In conclusion, we have developed a theoretical descrip- 
tion of the exciton-acoustic phonon scattering in disor- 
dered QWs which goes beyond the Markov approxima- 
tion. We have considered deformation potential cou- 
pling with acoustic phonons and realistic exciton states 
resulting from a microscopic simulation. A second or- 
der self energy approach is used for computing the ab- 
sorption spectrum from the exciton propagator. At 
the Markov level the absorption is a superposition of 
Lorentzian peaks with widths defined by Fermi's golden 
rule, which results from phonon-assisted scattering be- 
tween different states. At the full non-Markovian level, 
virtual scattering events on the same state are consid- 
ered as well (pure dephasing). These processes mani- 
fest itself in broad absorption bands around Lorentzian 
peaks (ZPL) whose widths are smaller than according 
to Fermi's golden rule. We have also discussed a semi- 
Markov approximation in which a diagonal self energy 
is employed, but the important frequency dependence of 
the diagonal term is preserved. This allows to reduce the 
numerical task from 0(N 3 ) to 0(N 2 ) 1 where N is the 
number of exciton states, while recovering all the impor- 
tant features of the complete solution. We have shown 
that the non-Markovian effect increases with a shorter 
localization length of the exciton wave function, and for 
low spectral density of states. A near-field setup could 
demonstrate these predictions experimentally. 
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APPENDIX A: INTERACTION MATRIX 
ELEMENTS AND COUPLING FUNCTION 

The matrix element for optical interband transitions 
between an exciton wave function ^f a (r e rh) and a light 
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field with momentum q and field envelope E(r) is given 

by 



dre- iqr £;(r)* Q (r,r) 



(Al) 



with the interband dipole moment /i cv of the band edge 
states. The factorization Eq. (1) and transmission normal 
to the QW plane (qii = 0) results in 



fi cv cj) ls (0)O eh J dRE(r) V» a (R) , (A2) 



with 0i s (0) as the relative exciton wave function at the 
origin, and O e h the confinement overlap integral O e h — 
J dz u e (z) u h (z). 

Within the same factorization Ansatz, the deformation 
potential matrix clement of Ref. 24 takes the form 



l af3 — 



hwc. 



2u 2 s p M V 



Oa#)q|, X 



D c K e (q z )x(q.\\/Ve) - D v K h (q z )x(<l\\/Vh) 

(A3) 

with deformation potential constants D a , sound velocity 
us, mass density ps, and normalization volume V. Fur- 
ther, we have defined Fourier transforms of the squared 
confinement and relative wave functions, and the COM 
overlap 

K a (q z ) = JdzulWe-*" 

X (qn) = JdptfMe-^rP (A4) 

(V^q,, = JdR^ a (H)e-^^n). 

In particular, for Gauss confinement states and a 
hydrogen-like relative motion we obtain: 



CX P ( -^(lzL a ) 2 



K a {q z ) 

X(q||) = (l + ( 9 ||a B /2) 2 )- 3/2 



(A5) 



where L e = 1.69 nm and Lh = 1.34 nm. These two values 
result from a variational determination of the confine- 
ment levels in the 5 nm wide GaAs/AlGaAs QW. Other 
material constants for GaAs are listed in Tab. I. 

Since the exciton states ip a (R) arc numerically gen- 
erated in Cartesian coordinates, we proceed with fast 
Fourier transformation along q x , q y , q z and integrate ac- 
cordingly Eq. (10) as 

J ar,( E ) = J dq x dq y dq z tl fj t-* X 

5{\E\-hu s JqJ+qJ+ql.) (A6) 



The delta function is used to integrate over q z . Special 
care is taken to deal with the resulting square root singu- 
larities at the remaining integration limits. Each matrix 
element t contains a square-root of phonon energy, thus 
J@ V (E) starts at least with a power of E 3 . 



APPENDIX B: ZPL WIDTH AND WEIGHT 



In the semi-Markov approximation, the absorption is 
given by 



a(uj) = Im^ 



m c 



(Bl) 



For evaluating the ZPL width, we expand in Eq. (Bl) 
the self energy around the pole e a . The real part at the 
resonance (polaron shift) ReSf M (e Q ) is small and can be 
dropped, but its derivative has to be kept as 



— S a 



(B2) 



LO — € a — iO 



The imaginary part is taken as ImSf M (e Q — iO) = 7^/2. 
Putting this into Eq. (Bl) we end up with an absorption 
lineshape 



which has a reduced ZPL weight of 1/(1 + S a ) and a 
reduced width 7 f M = 7^/(1 + S a ). 

At high temperatures, the phonon occupation can be 
replaced by n(E) — > E/ksT. Together with the Gaussian 
shape of the diagonal coupling factor (see Fig. 3) we can 
evaluate the Huang-Rhys factor Eq. (14) in closed form 



S a ^ JdE^E 3 Acxp 



A n E 



± = h I_ 

E* A Q 



Thus, the ZPL weight is approximately given by 



1 



1 + S a A a + bT 



(B4) 



(B5) 



This approximation is getting better for less localized 
states, which can be seen in the top panel of Fig. 7. 
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